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Abstract 

We derive an analytical expression for the distribution of the fcth smallest Dirac eigenvalue 
in QCD with imaginary isospin chemical potential in the Dirac operator for arbitrary gauge field 
topology v. Because of its dependence on the pion decay constant through the chemical potential 
in the epsilon-regime of chiral perturbation theory this can be used for lattice determinations of 
that low-energy constant. On the technical side we use a chiral Random- Two Matrix Theory, where 
we express the kth. eigenvalue distribution through the joint probability of the ordered k smallest 
eigenvalues. The latter can be computed exactly for finite and infinite N, for which we derive 
generalisations of Dyson's integration Theorem and Sonine's identity. 



1 Introduction 



It is by now well known how the spontaneous breaking of chiral symmetry in QCD leads to remark- 
ably strong predictions for the spectral properties of the Dirac operator in that theory. Based first 
exclusively on the relation to the effective field theory for the associated Nambu-Goldstone bosons 
at fixed gauge field topology [lj, an intriguing relation to universal Random Matrix Theory (RMT) 
was also pointed out [2J. It has subsequently become clear how these two alternative formulations are 
related, and all n-point spectral correlation functions have been shown to be identical in these two 
formulations to leading order in a 1/L-expansion (where L = y 1 / 4 gives the extent of the space-time 
volume V) OH]. This holds then also for individual distributions of Dirac operator eigenvalues [5]. 

If one seeks sensitivity to the pion decay constant F w it turns out to be useful to consider the 
Dirac operator of quark doublets with isospin chemical potential /x. Based on the chiral Lagrangian 
formulation [6], it has been suggested to use a spectral 2-point function of the two associated Dirac 
operators with imaginary isospin chemical potential. The advantage of imaginary chemical potential 
lies in the fact that the corresponding Dirac operator retains its anti-hermiticity. To leading order, 
all results can be expressed in terms of the simple finite-volume scaling variable (l = fiF^y/V. In this 
way, F n can be extracted from fits that vary \i and/or V. There is also sensitivity to F n in other 
observables that couple to chemical potential [3 [U EJ QUI E] ■ 

The leading-order chiral Lagrangian computations of ref. [6] have been given a reformulation 
in terms of a chiral Random Two-Matrix Theory in ref. [12]. In this way, all spectral correlation 
functions associated with the two Dirac operators T>\ and T>2, with respective chemical potentials \i\ 
and fi2, have been computed analytically in [12] for both the quenched and the full theory with Nf 
light flavours. It also includes all spectral correlation functions where the imaginary isospin chemical 
potential only enters in the Dirac operator whose eigenvalues are being computed, while the gauge 
field configurations are obtained in the usual way at vanishing chemical potential. In analogy with 
what is being done when varying quark masses away from the value used for generating the gauge 
field configurations we call this "partial quenching" . 

In an earlier paper [13] it was shown how all probability distributions of individual Dirac operator 
eigenvalues can be computed by means of a series expansion in higher n-point spectral correlation 
functions. In reference [13J an explicit analytical formula was also given for the lowest non-zero 
eigenvalue distribution for arbitrary combinations of flavours N± and N2 of the Dirac operators T>\ 
and T>2, respectively, at gauge field topology v = 0, in an approach quite close to that of ref. [H]. Two 
obvious questions remained open that we will answer in this paper: how to extend this to non-zero 
topology v > 0, and how to compute the distribution of the second, third or general feth eigenvalue as 
these are known for zero chemical potential [15]. However, the path of ref. [13] is not very suitable in 
particular for the derivation of the distributions of higher eigenvalues in a compact analytical manner. 
Our setup will follow closely ref. [15] at vanishing chemical potential. We shall present here a new 
formalism that immediately allows for the analytical determination of the distributions of these higher 
eigenvalues. The extension of the first approach [13J to v > for the first eigenvalue is presented in 
appendix \C\ as an alternative formulation and analytical check to part of our new approach. The 
benefit of our new results should be two-fold: while expressions for higher topology v > allow for 
an independent determination of S and F n from different lattice configurations, the expressions for 
higher eigenvalues should allow for a better determination using the same configurations as for the 
first eigenvalue. 

How much of the program proposed in [6] to determine F n on the lattice has been realised in the 
meantime? Based on a preliminary account [16] of ref. [13] , the expansion of the first eigenvalue was 
first used in simulations in [T7]. However, the question remained how large the finite- volume corrections 
to the leading order (LO) e-expansion are, in which the chiral Lagrangian-RMT correspondence holds. 
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In a series of papers this question has been addressed and answered: in |18j the next to LO corrections 
(NLO) and in [19] next to next to LO (NNLO) corrections in the epsilon-expansion were computed. As 
a result of these computations at NLO all RMT expressions for arbitrary n-point density correlations 
functions (and thus for all individual eigenvalues too) remain valid. The infinite volume expressions 
simply get renormalised by finite-volume corrections, one only has to replace S and F n by £ e // an d 
F ne ff in the corresponding RMT expressions. Here the subscript "e//" for effective encodes the 
corrections that match those computed earlier in pOJ and |21|, lllj . respectively. Only at NNLO 
non-universal, non-RMT corrections appear. It was further noticed in [181 119j. that the size of the 
corrections at each order depends considerably on the lattice geometry, in particular when using 
asymmetric geometries. In order to keep the NNLO corrections small, in [22 j the authors used a 
specific optimised geometry where they could apply RMT predictions and effective couplings at NLO 
only, and they obtained realistic values for S and F n from partially quenched lattice data for small 
chemical potential. Technically speaking £ was determined there from the first eigenvalue distribution 
at vanishing chemical potential, and F n from the shift of the eigenvalues compared to zero chemical 
potential. We refer to [22] for a more detailed discussion of these fits. 

Motivated by these findings we have completed the computation for the feth Dirac eigenvalues for 
all v > in the RMT setting, in order to have a more complete mathematical toolbox at hand. 

Our paper is organised as follows. In the next section we briefly define the notation and remind 
the reader of the definition of chiral Random Two-Matrix Theory. We introduce a certain joint 
probability density and describe how it can be used to derive individual eigenvalue distributions. We 
give the explicit finite- N solution here in terms of new polynomials and a new sequence of matrix 
model kernels. We take the scaling limit relevant for QCD in section [H and write out explicitly and 
discuss the physically most important examples such as partially quenched Nf = 2 results. Section 
[5] contains our conclusions and a suggestion for a quite non-trivial but important extension of these 
results. Because some of the relevant technical details have been described in ref. [13], we have 
relegated many of the mathematical details in this paper to appendices. In addition, in Appendix 
A we describe an explicit construction of the polynomials needed to compute the first eigenvalue 
distribution in sectors of non-trivial gauge field topology if one alternatively uses the method of ref. 

2 Chiral Random Two-Matrix Theory 

Before turning to the relevant Random Two-Matrix Theory, we first briefly outline the set-up in the 
language of the gauge field theory. We are considering QCD at finite four-volume V, and we assume 
that chiral symmetry is spontaneously broken at infinite volume. We consider two Dirac operators 
Dl,2 with different imaginary baryon (quark) chemical potential /Ui,2, 

Dx4 n) = [#(A) + ^ l70 ]4 n) = iX^[ n) 

D*$> = [0(A) + w>b]l4 B) = i>P4 n) ■ (2-1) 

When fi = Hi = —[12 this is simply imaginary isospin chemical potential, but we can stay with the 
more general case. We thus consider N\ light quarks coupled to quark chemical potential fj,%, and N2 
light quarks coupled to quark chemical potential ^2- Let us first consider the conceptually simplest 
case where N\ + N2 = Nf, and later comment on the changes needed to deal with partial quenching. 

In the chiral Lagrangian framework the terms that depend on are easily written down on 
the basis of the usual correspondence with external vector sources. Going to the e-regime of chiral 
perturbation theory in sectors of fixed gauge field topology pQ, the leading term in the effective 
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partition function including imaginary fi\^% reads [HI [23] (see also [24] for QCD-like theories) 

zi Nf) =l dU ( detU y e kVF2Tr[U,B)[rf,B]+^VTr(MiU+MUl) 
JU(Nf) 

In ([22D the N f x N f matrix 

B = diag(fj, l l Nl , /jl 2 1n 2 ) ( 2 - 3 ) 
is made out of the chemical potentials, and the quark mass matrix is 

M = diag(mi,...,mAr / ) • (2.4) 

The partition function (|2.2p is a simple zero-dimensional group integral. The leading contribution to 
the effective low-energy field theory at finite volume V in the e-regime is thus well known. 

We consider now the limit in which V — > oo while rh = mHV and p, = fiF w \^V are kept fixed. 
In this limit, to LO in the e-expansion, the effective partition function of this theory and all the 
spectral correlation functions of its Dirac operator eigenvalues are completely equivalent to the chiral 
Random Two-Matrix Theory with imaginary chemical potential that was introduced in ref. [12]. 
The equivalence for the two-point function follows from [6], for all higher density correlations it was 
proven in [4]. Therefore, since we have proven [13] that the probability distribution of the A;th smallest 
eigenvalue can be computed in terms of this infinite sequence of spectral correlation functions, we are 
free to use the chiral Random Two-Matrix Theory when performing the actual analytical computation. 

As already mentioned it has been shown in [18] that also to NLO in the e-expansion the random 
matrix expressions [12] for density correlation functions remain valid, when replacing £ and F n by 
the renormalised constants £ e // an d Fneff that encode the finite- volume corrections. Only to NNLO 
non-universal corrections to the random matrix setting appear. 

The partition function of chiral Random Two-Matrix Theory is, up to an irrelevant normalisation 
factor, defined as 



zf f) = f d$d* e -^(* t *+* t *) Y[ det[2?i + m h ] f| det[P 2 + m h ) (2.5) 

J /l = l /2 = 1 

where T>\ 2 are given by 



v ^-\ ist + j^t J • (2 - 6) 



The operator remains anti-Hermitian because the chemical potentials are imaginary, as shown ex- 
plicitly. Both $ and $ are complex rectangular matrices of size iV x (N + v), where both N and v 
are integers. The index v corresponds to gauge field topology in the usual way. The aforementioned 
correspondence to chiral perturbation theory holds in the following microscopic large- ./V limit: 

lim zi Nf) = zi Nf) with rh = 2Nm , p, = V2N y, . (2.7) 

N— >oo 

In the framework of chiral Random Two-Matrix Theory it is particularly simple to consider the 
situation corresponding to what we call partial quenching. Here one simply considers eigenvalues of 
one of the matrices, say D\, that then does not enter into the actual integration measure of (|2.5p by 
setting N\ = 0. In the language of the chiral Lagrangian, this needs to be done in terms of graded 
groups or by means of the replica method. 
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Referring to ref. |12] for details, we immediately write down the corresponding representation in 
terms of eigenvalues xf and yf of T>\ and V 2 , respectively, 



Z " N/) = / X\dxid yi ri Nf \{x},{y}) , (2. 
Jo 7=1 



up to an irrelevant (mass dependent) normalization factor. The integrand is the joint probability 
distribution function (jpdf), which is central for what follows: 

N I Ni N 2 

ri Nf \{x}, {y}) = [J (x iyi y +1 e- N ^+^ ]J {x\ + m 2 ;i ) J] (»? + m D 

i=l \ /i=l h=l 

x A A r({x 2 })A 7V ({y 2 }) det [J„(2diVx i y i )] • (2-9) 

l<i,j<N 

Because the integration in eq. (|2.5p was over $ and ^ separately, the matrices now become coupled 
in the exponent. The corresponding unitary group integral leads to the determinant of modified I- 
Bessel functions, and removes one of the initially two Vandermonde determinants, which is defined as 
Aat({x 2 }) = r[j>j=i( x j ~~ X D- The precise connection between the constants and ^1,2 is given by 

ci = (1 + vh/S 2 , c 2 = (1 + 5 2 , 
d = (1 + mn 2 )/5 2 , 5 = [is -m, 
1-r = d 2 /(c lC2 ) , (2.10) 

where the latter is defined for later convenience. We need the joint probability distribution to be 
normalised to unity, which is done trivially by dividing by Z„ Nf (cf. eq. 12. 8p ). 



3 The kth Eigenvalue at Finite- TV for Arbitrary v > 

We now follow the derivation of ref. [15] rather closely. We are able to do that because we focus 
here on the distributions of individual ^-eigenvalues only - which are those we may partially quench. 
For that purpose it is convenient to first consider the joint probability distribution of the k smallest 
x-eigenvalues, ordered such that < x\ < x 2 < . . . < 



ni Nf \ Xl ,...,x k ) = • / dx k+1 --- dx N T\dy l vl Nf \{x},{y}). (3.1) 

Z { u Nf) (N -k)\Jx k Jx k Jo f = \ 

This quantity is then used to generate the kth x-eigenvalue distribution through the following inte- 
gratiorQ 

pf s ' u \x k )=[ dxi [ dx 2 ... [ dxk-x Sli Nf \xi,...,Xk) . (3.2) 

JO J XI J x k _ 2 

Note that for k = 1 no integration is needed, and Pi f ' (x±) = fi£ ■ (si). 

The computation of mixed or conditional individual eigenvalue distributions, e.g. to find the joint 
distribution of the first x- and first y-eigenvalue, remains an open problem. 



1 Compared to [15] we are already working with squared variables here. Translating to that picture the integration 
bounds in eqs. (|3.1[) and (|3.2|l remain the same as in [15] . 
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We next proceed as in ref . [13] , and integrate out all y-eigenvalues exactly. Because of this we note 
that in eq. ()3.1|) we can replace the determinant over the Bessel functions by AM times its diagonal 
part, after having made use of the antisymmetry property of A]y(y 2 ). After inserting a representation 
of the Bessel function in terms of a factorised infinite sum over Laguerre polynomials (see eq. (B.7) 
in [12]), we get 



f oo N fOQ N I Ni N 2 

/ n fa ^(m. {y}) = w n w n & + o n ^ + 

Jo i=i Jo i=i V / 1= i / 2 =i 



m%) A N ({x 2 }) 



N 



^iv({ y 2 })II ( W^^W^e-^^-W) £ + ffi K {Nr Cl x 2 W n% (Nrc 2 y 2 )) , 



i=i 



rij=0 



(3.3) 



where the Laguerre polynomials L v -{NTC 2 y 2 ) now appear with their corresponding weight function 
y2i/+i e ~NTc 2 y 2 j ue ^Yi e identity used. Next we include the set of N 2 masses, {1712}, into Ajy({y 2 }) 
to form a larger Vandermonde determinant of size N + N 2 , and then replace it by a determinant of in 
general arbitrary Laguerre polynomials normalised to be monic 



N N 2 

A N ({y 2 })HH(y 2 + m% 

i=i / 2 =i 



LUNrc 2 (im f2=1 ) 2 ) 

L D Q (NTC 2 (im N2 ) 2 ) 
Ll{Nrc 2 y 2 ) 

Ll{Nrc 2 y 2 N ) 



(Ntc 2 



) l+ N2 -i L p N+N2 _ 1 (NTc 2 (im f2=1 ) 2 ) 



{NtC2) N+N 2 -1 L% + N 2 -l ( Ntc 2 ( im N 2 f 
(Nrc 2 )^ J + N 2- 1 L N+N 2 -l( NTC ^yi) 

JW^ 2 yr+N 2 ^L u N+N2 _ l (NTc 2 y 2 N ) 



A N2 ({(im 2 ) 2 }) 



(3.4) 

Here the index of the Laguerre polynomials v is arbitrary. The monic Laguerre polynomials relate to 
ordinary Laguerre polynomials in a ^-independent manner: 



L»(x) = (-l)"n! L» n {x) = £(-1) 



n+j _ 



n!(n + v)\ 



j=0 



{n - j)\{v + 



n-l\ 



(3.5) 



In eq. (|3.4p the inverse powers {Nrc 2 y~ l can be taken out of the determinant. Inserting this back into 
eq. (13. 3f) for v = v we can use the orthogonality of the Laguerre polynomials in the integrated variables 
2/1 ... N j killing the infinite sums from the expanded Bessel functions. The Laguerre polynomials in x 
thus replace those in y inside the determinant, times the norm from the integration. We obtain 



poo N 

/ Y[d yi vl Nf \{x},{y}) = 
Jo i=i 

M(Nd) N,/ T N ( u+ V U?Jo N *~ X Q- ~ T ) j (NTc 2 y 
A N2 ({(im 2 ) 2 }) 2 N (Nrc 2 ) N ^+ 1 ) 



N 

n 

i=i 



x 2iy+l e -NT Cl xf 



H(x 2 + m 2 fl )\A N ({x 2 }) 
/i=i 
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L v {NTC 2 {im h=l ) 2 ) 

L»{NTc 2 {im N2 ) 2 ) 
L^NTaxj) 



L v q (Ntc iX 



N) 



{1 - T) N+N 2 -i L p N+N2 _ l (NTC 2 {im h=l ) 2 ) 

{1 _ T) N + N, 2 -i L v N+N2 ^ l {NTc 2 {im N2 ) 2 ) 
Ln+n 2 -i(Nt Ci x 2 ) 

L^^iNrcxxjy) 



(3.6) 



after taking out common factors of the determinant. The determinant in eq. (I3.6p . which we call 
-Dtv_i_tv 2 , can almost be mapped to a Vandermonde determinant, using an identity proved in appendix 
A in \m 



DN+N 2 ({m 2 };{x 2 }) 



L v o{W 2 N2 ) 
1 

1 



T N+N 2 -1 r a 2 

( 1 _ t )JV+jv 2 -i ^N+N 2 -l V T 1V1 f 2 =l I 

T N+N 2 -1 ~ J 2 

(1_ t )N+n 2 -i-^n+N 2 -1\ t 1V1 N 2 ) 



X 



2(N+N 2 -1) 
1 



X 



2(N+N 2 -1) 
N 



(3.7) 



where we have defined 



Mj 3 = NTc 2 (im f2 ) 2 and X 2 = Ntc x x 2 . 



(3.8) 



This fact can be used below to perform the N — k remaining integrations in the generating quantity 
ilje, after inserting eq. (|3.Tj) into eqs. (|3.6p and (|3.ip . This leads to 



j>i>i 



2 2a 
, 



/ • • • / dx N yi ~ x ?) n 

•^ Xfc j>i>k+l j=k+li=l 



2 2a 



i=l 



/l=l 



(3.9) 



where we have split the Vandermonde determinant An({x 2 }) into integrated and unintegrated vari- 
ables, and defined the following constant 



C = 



(N\) 2 {Nd) Nv T N ^ n™ 2 "^ 1 - tY{Ntc 2 )-i 
zl Nf) (N-k)\ 2 N (Nrc 2 ) N (^A N2 ({(im 2 ) 2 }) 



(3.10) 



We can now change variables Xj 



x 2 . for j = k + 1, . . . , N, and then perform the shift 



to obtain integrations dzj in eq. (|3.9|) : 



nl N f\ Xl ,..., Xk ) = c ]J(x 2 -x 2 )ll 



x 2v+l e -NT C1 x 2 t 



-N(N-k)Tdxl 



j>i>l 



L n [** 



N1 \ 1 

\{(x 2 + m 2 h ) 
' - J \ /i=i / 

fc-l iVi 
e ""«** ( z . + 4)^ JJ(^ + x | _ x 2 } JJ + x 2 + m 2 ;i 

t=l / x =l 



V 



X 1 I (^j -^0 ^A r +A r 2 ({ ?7i 2}; 2; Ii---> :l; l) 2; fe+l+ 2; L---^A r + :E l) • ( 3 - n ) 
j>i>fc+l 
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1 extra mass terms of flavour-type "1", in addition to the N\ 

w(z) 



We thus obtain an integral with v + k 
shifted masses. The weight 

l e~ NTClZ (3.12) 

is now of Laguerre-type corresponding to a fixed topological charge of D = 1, irrespective of the actual 
topological charge v of the given gauge field sector we started with. We will therefore call v spurious 
topology. Compared with the corresponding derivation in case of vanishing chemical potential [15], 
this can be seen to differ by one unit, compared to spurious topology v = 2 at vanishing in [T5] . 
The reason for this difference is easily traced to the different integration measure for the x-eigenvalues, 
which has one power less in the Vandermonde determinant compared to the case of vanishing chemical 
potential. It is an interesting and quite non-trivial check on our present calculation that we recover 
the results of reference [15] in the limit of vanishing chemical potential. In particular, the shift from 
spurious topology v = 1 to spurious topological charge v = 2 in the integration measure will now arise 
due to recurrence relations of Laguerre polynomials. Some details of this will be given below. 

When replacing the Vandermonde determinant in the variables Zj as well as Dn + n 2 by a determi- 
nant containing Laguerre polynomials we thus choose polynomials L^(Ntciz) in order to be able to 
exploit the orthogonality properties with respect to the measure w{z) eq. (I3.12p . 

For the new masses times Ajv-fc({^}) this is an easy task. We can include them into a bigger 
determinant of size N — k + N% + v + k — 1, following the identity eq. (13. jj) . Here we replace the N 
variables yf by N — k variables Zi, and the set of N2 masses by the following set of N% + v + k — 1 
masses: 



and likewise we define 



m 



m' 2 

m h 
12 

m Nl+j 
12 



m 



+ 4 

x 2 



for /1 = 
for j = 1 
for i = l. 



l,...,iVi 
, . . . , k — 



1 



(3.13) 



M 



12 



NTCi(im'j) 2 for j = 1, 



, JVi + v + k - 1 



(3.14) 



For computational simplicity we first set the v degenerate masses to be different by adding small 
pairwise different constants, e|, and then set ej = at the end of the computation. Also we may chose 
spurious topology v = 1 in eq. fl3.4f) . The prefactors in front of the Laguerre polynomials inside the 
determinant can be taken out. 

To express the determinant Dn+n 2 °f the shifted arguments in eq. (I3.11j) in terms of Laguerre 
polynomials requires a bit more algebra: 

D N+Na ({ml};xl,.. 



L V ^M 2 ) 



1 X}~i Zfc- 



1 + x \i 



Z N 



x k) 



2^=o (i_ T )(^ 



E 



N+N 2 - 




-1 t 1 furl 
{XI 



M 



n 



,JV+JV 2 -1 



x 2 k) N+N 2 -l 



Z fc+1 



z 



N+N 2 -l 
N 



2 Of course, the additional pieces due to the ^-integrations are what ensures equivalence to those corresponding 
one-matrix model results in the limit /ij — > 0. 
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;(m 2 



Q u o(M 2 



ya=i" 

iv 2 ) 



LI(Zn) 



q u N+N2 -x(M 2 N2 ) 

±J N+N 2 -l 



^iV+JV 2 -l(^fc+l) 



f 1 

N+N 2 - 



(Zn) 



(3.15) 



where for convenience we have defined M'jf +u+k 



0, as well as 
Z k = Nrciz k . 



(3.16) 



In the first step in eq. (|3.15p we have used the invariance of the determinant to undo the shift in xi 
of the Zj variables. This leads to a shift in variables Xi and to linear combinations of the Laguerre 
polynomials in the N2 masses. In the second step we have added columns from the left to the right 
to replace monic powers in Zj and M' 2 by polynomials L\. Because the determinant D^+n 2 is n °t 
an invariant Vandermonde this leads to a further sum in the first N2 rows, invoking the following new 
polynomials 



1 

1=0 \ ' 



(3.17) 



The form given here is derived in Appendix [A] using identities for Laguerre polynomials. For later 
purpose, we note already that in the limit of zero chemical potential, i.e. in the limit r — > 0, we obtain 
Laguerre polynomials of shifted mass from the 



limC(M 2 2 ; 

T— W 



) n n\ 



Y,m-Nml)L 

j=0 



v 
n-j 



(-Nxt) = Li(-N(ml + xl)) 



(3.18) 



In this way we recover, after the use of a few identities for Laguerre polynomials, the results of ref. 
|15j in the limit of vanishing chemical potential. 

We now proceed with the integration over the variables z^+i, ■ ■ ■ , zjy in eq. (|3.1ip . Using the 
rewriting discussed above, we have: 



(AM 



A'l 



X\, . . . , X/jJ - 

n 

j=k+l 



c n 'j m 

j>i>l i=l 

N+Ni+v-2 

y ■ p-NrciZj 



) n ( ntc ^ 
3=0 



r2\ 



Ll(M[ 



f 1 



>2\ 



L N+N!+u-2( Z k+l) 



\{M 2 



Ll{M' 2 



N 2 ) 



Ll{M' 2 i+v+k ) 
L (Z k+ i) 

l o( z n) 



2~{N-k) e -N(N-k)rcixl 
Ajv 1+I)+t: -i({(im') 2 }) 



<&+N 2 -l{M} a=1 ) 



qN+N 2 ~i 

1 

J N+N 2 - 



-i(M ,2 1+v+k ) 

L N+N 2 -l( Z k+l) 



T 1 

^N+N 2 



l n+n 2 -i( z n) 
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(3.19) 



This expression can be simplified somewhat by noting the identity 



k-Nt+v+k-iiiiim') 2 }) 



A JVl ({(im 1 )2})A,({(z e )2})nf =1 m^ 



'2v 



;i + 0(e)) . (3.20) 



The Vandermonde determinant Ajv x ({(imi) 2 }) also occurs as a factor in the partition function, and 
thus cancels. 

The integrals over Zj,j = k + 1, . . . , N can be performed, again exploiting orthogonality properties 
of the Laguerre polynomials L^z) with respect to our weight function f|3. 12|) . This is done by means 
of a generalisation of the original Dyson Theorem (see e.g. [25]), now for two determinants of different 
size, of which different sets of entries are not even of Laguerre-type, and not integrated over either. 
The needed generalisation of this theorem was essentially provided in ref. [26], and we only need a 
slight extension of this more general theorem here. We have relegated the proof of the new theorem 
to Appendix [Bj and will only quote the result here. 

To express the result, we define two kernels (both corresponding to the kernel K of Appendix IB1 
for different values of the indices) 

K{M[\M?) = ^-1^L}(M( 2 )L}(M^ 2 ) , K»(M> 2 ,M 2 ) = £ J^L)(M' 2 )^(M 2 ) , (3.21) 

with A = min(A + N\ + v — 2, N + N 2 — 1), and the associated matrix of size (N 2 + k) x (JVj. + u + k — 1) 

/ K»{M[\M% =1 ) ••• K»(M> N \ +v+k _ v M% =1 ) \ 



K 



K»(M[ 2 ,M 2 N2 ) 
K{M' 2 ,M' N \ +V+1 ) 

\K(M[ 2 ,M' N 2 i+u+k ) 



K v {M' 2 i+u+k _ v M 2 m ) 
K(Mg 



,2 M ,2 



+v+k—V lv± Ni+v+l> 



(3.22) 



K(M' 2 i+u+k _^M' 2 i+u+k )j 



We can now give the answer for Q u } , where we have to distinguish 3 different cases, depending 
on the respective sizes of the two determinants in (|3.19|) . 

i) The simplest case is N 2 = N\ + v — 1 (this implies that tC is square). Applying Appendix IB1 to 
the integral of (|3,19p and using eq. (|3.20p we obtain to leading order in mass difference ey. 



„2u+l rrfc-l 



.X k ) 



c- 



-Nrcixh 2 



Xt - X; 



2^^A Wl ({(fm 1 )2})A,({(z6) 2 })nX 1 =1 ^r 



{N TCl ] 



-2(N-k) 



N+N ± +u-2 



x e -m-"+i)rc lX l ] [ (AT TCl )-i J] / (j + l)\j\\(N - k)\(-)^ +N2+ ^ (N -^ detIC 
3=0 j=(A ' 

(3.23) 



A typical example for this case would be N- 



A, 



17=1. 
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ii) For N 2 > N\ + v — 1 the expression is the same, but with additional columns in the determinant 
in the second line of eq. ()3.23[) : 



^N+Ni+u-l(^Ni+v+j) "' ^N+N 2 ~l^'N 1 +v+j) 



, f 2 = l,...,N 2 , j = l,...,k 



(3.24) 

This case typically occurs for partial quenching N\ = 0, N 2 = 2 and small topology v = 0, 1, 2 

Hi) Similarly, for N 2 < N\ + u — 1, which typically occurs for higher topology, the determinant in 
eq. (|3.23|) is replaced by 

IC T L^ +iV2 (Mj 2 ) ••• L^ +JVl+1 ,_ 2 (Mf) |, j = l,...,N 1+ u + k-l. (3.25) 

We note that all three cases contain degeneracies for v > 1, and that a Taylor expansion must be 
performed when sending €j — > 0. 

Finally we are ready to compute the feth individual eigenvalue probability distribution (x) by 
means of inserting the above results into eq. (I3.2p 

Pk Nf ' U \ x k) = J dxi dx 2 --- dx k -i£ll Nf \xi, . . . ,x k ) . 

An alternative, equivalent formulation for the first eigenvalue Pu}-i{x) at arbitrary u follows from the 
derivative of its cumulative distribution given in appendix [Cj where we generalise the approach of [13J 
to v > 0. 

We refrain from giving further explicit examples for finite iV and instead turn to the large- iV limit. 

4 The Microscopic Scaling Limit 

Having obtained the explicit solution for any finite N, we are now ready to take the appropriate 
microscopic iV — > 00 scaling limit. In the language of QCD this corresponds to a finite- volume scaling 
in V, where V is the space-time volume. It is only in this limit that we expect to obtain universal 
results which do not depend on having chosen Gaussian measures for the two original matrices $ and 
fy. We remind the reader that the appearance of Laguerre polynomials in the finite- N solution is 
directly linked to these two original measures having been chosen Gaussian. As we have seen, there 
are numerous instances where our derivation makes explicit reference to specific identities for Laguerre 
polynomials, and also the measure factor was specific to these Gaussian integrals. It is therefore a 
quite non-trivial task, and an interesting challenge, to generalise the universality proof of the last 
reference in [2] to this more general setting. The universality has been implicitly checked by the 
equivalence proof of [3], but an explicit proof of universality directly in the framework of the chiral 
Random Two-Matrix Theory remains to be constructed. 

In QCD-terminology, we keep rhi = rriiT^V and fi = fiF^y/V fixed, while we take V — > 00. In the 
language of our chiral Random Two-Matrix Theory, we take the — > 00 limit while keeping 

obi = 2Nxi , rhi = 2Nmi , fi\ 2 = 2A>f >2 (4.1) 

fixed (we also scale the masses €j in the same way, before taking them to zero at the end of the 
calculation). In addition, we introduce the following relevant quantity, the difference in rescaled 
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chemical potential, 5 = fa — Al- We will follow |13j closely, without giving a detailed derivation. 
Defining the following continuum indices 



t = j/N, and r = l/j 



(4.2) 



we replace the sum by an integral, ^2f =0 N 1 — > J dt, and correspondingly for index I. In addition, 
we use the identities 



lim 77 

N,j,l^oc (1 - t) 1 



exp 



-rtd 2 



lim j~ v L\{M = -Nram 1 

N,j,l^oc 



4r 



m 



v/2 



I u {Vrtm) . (4.3) 



The scaling limit of the kernel K is obtained using the Christoffel-Darboux identity: 

Ksirti^m'z) = lim -Lr n (M[ 2 , M' 2 2 ) 

1 Lj v+1 (M{ 2 )Lj v (M^ 2 ) - L\, +l (M' 2 2 )L),(M[ 2 ) 



N™oc M' 2 -M' 2 2 
8 m! x Iq (rhi ) I\ (m 2 ) — m 2 Iq (m 2 ) I\ (m^ ) 



1"*2 



* 1 2 " 1 1 

— m 2 



(4.4) 



The scaling limit of the new polynomials q v N requires a little more care. For u > one needs to 
treat the terms corresponding to indices I = N — z/+ 1, . . . , iV in (|3. 1T[) separately, due to the presence 
of Laguerre polynomials with negative index^l [13]. For these "anomalous" terms, the scaling limit is 



*)= *a i" 1 £ Tf—Ty L r ( M2 ) L 7-z y 

N.i^roo * — ' ( 1 — T) 



lim \ 



-l/~p-l 



7^ 



-i 



p=0 



P 3-P 



E (-r' / "- p - 1 



E 

p=0 



1 ftxlY d u - p - 



Q r V-p-l 



r=l 



4r W 2 



j-l-pj (1 -t) 



1 



^(v rim) 



(4.5) 



Here we write out and use the fact that the binomial weight kills the terms of order lower than 

v — p — 1 of the expansion in r. Note that naively q u SA is proportional to j" , but that the pattern 
of cancelation is exactly such that the limit is finite. The notation in eq. (|4.5|) is chosen such that for 
v = the sum is void, and thus the anomalous term is absent for v = 0, q v £^ = 0. 
For the remainder of the terms, the sum turns into an integral: 



<lSR\™> 



dr ea 



rt8 2 



1 - r 



v/2 



I v (y/rim)I v (y/(l - r)tx) . 



The final scaling limit of the new polynomials it then the sum 

q» s (m; 1) = lim tl q »( M 2 ) = q» R (m; t) + q v SA {rh; t) 
iVj->oo f.J J 



(4.6) 



(4.7) 



For integer v > and fixed j, z, the limit lim n ^oo n"L v (z/n) is only finite for j > v. 
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Let us give the limiting polynomials for the first few topologies v = 0, 1, 2 as examples: 

qT\rh-t) = ^°(m;t)= / 1 dre^ 2 7 (v^m)/ (7(1^7jtx) , (4.8) 



o 

q» =1 (m;t) = ^ [ % dr e\ rtP ( h(Vr~tm)h{J (I - rjix) + e^* -^h(y/im), (4.9) 



q u s = 2 (rh;t) = [ t dr -!— I 2 (V7trh)I 2 {^ {1 - rjix) 

m/ io 1 - r 



+ \^±^I 2 (Vtm) + -l-I^Vim^ . (4.10) 

For comparison see eq. (|C.12p for the limiting polynomials of the previous approach |13| for the first 
eigenvalue extended to v > 0. Because of technical reasons the anomalous terms appear there already 
at v = 0. 

The last building block we need is the new kernel K U (M' 2 ,M 2 ), eq. (f!T2"Tj) . It contains both the 
regular and anomalous terms q u s = q$ R + q$ A inside the integral (originating from the sum), but no 
further complications occur. Hence we have the scaling limit 

K%(m',m) = lim -L^(M' 2 , M 2 ) = A f dt Vth(Vim')q u s (m;t) . (4.11) 
N— >oo N 771 Jq 

Using these building blocks we can calculate the main quantity of interest, the scaling limit of the 
eigenvalue distributions, defined as 

Psk' U (^k)= / dxi dx 2 ... dx k _i n ( s J'(xi,...,x k ) , (4.12) 

JO J x\ J Xk-2 

...M = Jim (2N)-^ (xt = |L . . . , Xk = *£j . (4.13) 



with 



4.1 Example partial quenching 

As pointed out in [5], there is no 5 dependence in the fully quenched (N\ = N 2 = 0) theory. The 
simplest case to consider is then also what is probably the physically most relevant situation, namely 
the partially quenched case, where N\ = 0. On the lattice gauge theory side, this corresponds to 
generating the gauge configurations using N 2 dynamical quarks with zero chemical potential, and 
then looking at the spectrum of the Dirac operator T>\ with jjL\ = —5. 

In the simplest case N 2 = 1, v = and k = 1 (case ii), the matrix K, is empty, and (|3,24p is just 
a 2 x 2 determinant of polynomials, with M[ = 0. The two columns are degenerate, so we perform a 
Taylor-expansion in t. With the abbreviation rh = rhf 2= i, the limiting distribution is then found to 
be 



(o+i,o)/. \ e 2 r _i~2 



i«5 2 

(4-14) 



'ol 

-i.5 2 - 1 2 



q° s (rh;t = l) d t q° s (m;t) 
1 



t=i 
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in agreement (after partial integration) with the expression in [§J. By introducing primed masses, and 
changing the index of the (scaling limit of the) new polynomials qs and the partition function, we 
obtain the higher topology distributions which are new. 

For v = 1 (case ii), we can set e\ = in eq. (|3. 13|) from the beginning, and we find 



mx\e 



!*2_ 
" 4 1 



¥ 2 



16Ji 



m) 



1 



K s (x!,0) 



(4.16) 



where we note that the singularities of K(rh[, rh' 2 ) at rh 2 — > rh' x and rh 2 — > are removable. Increasing 
the topology further the primed masses become degenerate, so we take the scaling limit with the 
(scaled) e^-'s to be finite, and then let ej — > 0. The Vandermonde A„({(ie) 2 }) of (|3.23|) ensures that 
this limit is nontrivial. We then find for v = 2 (case i), 



(0+1,2),. ^ 

Psi 



* 2 -4 

m x\e 



4 1 



¥ 2 



64J 2 (m) 



Ks(xi, rh) drh'k 2 s {xi,m) 
K s (x u 0) d^KsixuO) 



(4.17) 



and for v = 3 (case Hi) 



m x\e 



4^1 



l^ 2 



64/ 3 (m) 



-RT|(fi, m) 
if 5(^1,0) 



^1' 



^1, 



9, 



'-^1(^1 



dl'Kl(xi,rh) 
dl,K s (x u 0) 



(4.18) 



where by dm' we mean the derivative with respect to the first argument (which, incidentally, is not the 
same as c^). Our new results for v = 1, 2 and 3 are illustrated in figure [TJ For alternative expressions 
in an equivalent formulation see appendix O The first eigenvalue distribution alone is in general not 
sufficient to fit both low-energy constants (LEC) S and F^, see the discussion in [22]. The benefit 
of computing the first eigenvalue distribution at higher topology here is to be able to fit both LEC 
independently from different lattice configurations, where typically v = 1 offers better statistics than 
v = 0. 




Figure 1: The first eigenvalue in the partially quenched case with N± = and ./V2 = 1 at mass rh = 3, 
shown for 3 different topologies and 2 different values of 5 = 1 (left figure) and 5 = 3 (right figure) 
The left, middle and right curves correspond to topology v = 1, 2 and 3, respectively. The effect of 
the increasing number of exact zero-eigenvalues can be seen nicely. We also show the corresponding 
spectral density from [12] for comparison (see eq. ([C.18P in appendix [C]). The first eigenvalue always 
follows the density almost up to its first maximum for all parameter values, which is an important 
consistency check. 
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Continuing with the new expressions for the kth lowest eigenvalues for Nf = + 1 with k > 1 , we 
will, for simplicity, focus on v = 0. Results at higher topology are easily obtained following the same 
approach as above. For k = 2 and k = 3 the joint probability distribution is found to be respectively 







(o+i) 

50 



(xi,x 2 ) 



e 2° 
16/ (m) 



x\X2(x2 — xf)e 



Ks(x 2 i,x 2 i 
K s (x 21 ,Q) 



q° s (m;t 
1 



1) 



d t q%{m\t) 

h(x 2 i) 




t=i 



. (4.19) 



and 



^SQ i.xi,x 2 ,x 3 ) 



e 2° 
1 28 J (m 



-x\X2Xz\x z — x 1 )[x 3 — x 2 )e 



4*3 



K° s (x 31 ,m) K° s {x 32 ) q° s (m-t = l) d t q%(m;t)\ t _ 

Ks(x 31 ,x 31 ) K s (x 32 ,x3i) ^h(x 31 ) h(x3i) 

Ks(x3i,x 32 ) K s (x 32 ,x 32 ) j^h(x 32 ) h{x 3 2) 

K s {x 31 ,0) K s (£32,0) " 1 




Figure 2: The first, second and third eigenvalue k = 1, 2, 3 in the partially quenched case with JVi = 
and N 2 = 2 at masses rh\ = 3 and m\ = 4, all at topology u = 0. We display 2 different values of 
5 = 1 (left figure) and 5 = 3 (right figure) to illustrate the influence of chemical potential. The left, 
middle and right curves correspond to k = 1,2 and 3, respectively. Also shown is the corresponding 
spectral density from [12] (see eq. (1C.24I) in appendix [C|) . It is compared to the sum of the first three 
eigenvalue distributions, following the density identically up to the third eigenvalue for all parameter 
values. 

The extension of these results to higher values of N 2 is straightforward. We will restrict ourselves 
here to N2 = 2 and v = 0. Defining the normalisation constant 



7^(0+2) 



Jo (mi) mill (mi) 
io(m 2 ) m 2 ii(m 2 ) 



(4.21) 



we have 



x\e 



lx.2 

"41 



AA°+2) 



d t q° s {mi;t) 
d t q° s (rh2;t) 



t=i 



d?q° s (m 2 ;t) 



t=i 
t=i 



(4.22) 
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Note that an alternative representation for the integral of this quantity, it cumulative distribution, 
was given in eq. (4.44) in |13j . 

For k = 2 and k = 3 we obtain the following new expressions 



n 



(0+2) 
50 



(X!,X 2 ) = - 



xix 2 m 



xi)e 4 2 



8AA(°+ 2 ) 

K° s (x2\,rhi) 
x Ks(x2i,m 2 ) 

Ks(x 2 i,x 21 ) 
K s (x 21 ,0) 



q° s {rh 2 ;t = 
1 



and 



n 



(0+2) 
50 



(xi,x 2 ,x 3 ) 



xix 2 x 2 (x1 — x\){x\ — x 2 )t 



-ir 2 
4 X 3 



K^(x 31 ,m 2 ) 
Ks{x 31 ,x 31 ) 
Ks(x 3 i,x 32 ) 



26^(0+2) 

K%(x 32 ,rhi) 
K%(x 32 ,m 2 ) 
Ks{x 32 ,x 3 i) 
Ks(x 32 ,x 32 ) 
K s (x 32 ,0) 



q%{mr,t 
q° s (m 2 ;t 



■r-M 
2 

■ r-,',2 



h{x 3 i) 
h{x 32 ) 



d t q° s {rhi;t)\ t=1 

kqlim^t) 1 1=1 

h{x 21 ) 




a 



d t q Q s {mi\t)\ t _ 
d t q° s (m 2 ;t)\ t __ 
h(x 3 i) 
h{x 32 ) 




Q2„0 



%q° s (m 2 -t)\ 




t=i 
t=i 



d^q%{fhi;t)\ t _ 
dh° s (^2;t)\ t _ 
^h(x 31 ) 
?fh(x 32 ) 




(4.23) 



(4.24) 



after inserting them into eq. (|4.12p . We illustrate the distributions of individual eigenvalues following 
from these equations in figure [2] for non-degenerate masses. In the case of equal masses, two rows 
become degenerate in both the numerator and denominator, necessitating a Taylor expansion. Figure 
[2] offers a further graphical consistency check for our results, when comparing the sum of the individual 
eigenvalues to the actual density. The two curves nicely agree almost up to the third local maximum 
for all parameter values. 



4.2 Two light flavours 

For Aq = N 2 = 1, v = the distribution of the first eigenvalue is given by 



(1+1,0) \ 

Psi ( x v 



x\(rh\ + xf)e' 



-if 2 

4 1 



8AA1+ 1 ) 
■i 



K^(m', rh 2 ) q° s (rh 2 ; t = 1] 
K s {m',0) 1 



A^ 1+1 ) = / die* 15 I (Vimi)I (Vim 2 ) , 
J o 



(4.25) 



(4.26) 



with rh\ = rhf 1= i, rh 2 = rhf 1=2 and m! = ymf + x 2 k . For its cumulative distribution see eq. (4.33) in 

The expressions for distribution of the second and third eigenvalue which are new results follow in 
a similar fashion from 



9. 



(i+i) 
50 



1*2 



(xi,X 2 ) 



x\x 2 (x2 - x\){rh{ + x A 2 )e 4^2 
6W( 1+1 ) 



K° s (m' ,m 2 ) K Q s (x 2 i,m 2 ) q° s (rh 2 t = I) 
K s (m',x 2 i) K s (x 21 ,x 21 ) ^h(x 21 ) 
K s (m',0) K s (x 2U 0) ' 1 



, (4.27) 
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and 



n (l+i)/„ , * X 
il y so '(xi,x 2 ,x 3 ) 



^1X2X3(^3 — xf)(x^ — x\)(rh\ + £3) 



e 43 



29^(1+1) 

K^(m',m 2 ) K^(x 31 ,m 2 ) k a s (x 32l rn 2 ) q° s (m 2 ;t = 1) 

K s (m',x 3 i) K s (x 3 i,x 31 ) K s (x 32 ,x 31 ) 

Kg(m', X32) ^5(^31,^32) ^5(^32,2:32) 

tf s (m',0) K s (x 3 i,0) i^(x 3 2,0) 



^/l(x 32 ) 



. (4.28) 



The corresponding figures are shown in figure [3j 




Figure 3: The first, second and third eigenvalue k = 1,2,3 for two light flavours with N± = 1 at mass 
rh\ = 3 and N 2 = 1 at mass rh\ = 4, all at topology v = 0. We display 2 different values of b = 1 
(left figure) and 5 = 3 (right figure) to illustrate the influence of chemical potential. The left, middle 
and right curves correspond to k = 1,2 and 3, respectively. The corresponding spectral density is 
taken from from [12] (see eq. (|C.22|) in appendix [Cj. Again we find an excellent agreement between 
the density and the sum of the three individual eigenvalues for all parameter values chosen. 

Both figures [2] and [3] illustrate the influence of chemical potential and thus the possibility to 
determine F n from individual eigenvalue distributions. In particular the shape of the first eigenvalue, 
but also the positions of the maxima of the second and third eigenvalue change considerably when 
comparing the left and right plots, in particular in the partially quenched setting in figure [21 Let us 
briefly compare to how £ and F n were determined from lattice data in ref. [22]. In order to avoid 
two-parameter fits, there S was first determined at fj,± 2 = from the first eigenvalue alone. Then, 
the two-point density correlation function p(x, y) of finding one x- and one y-eigenvalue was expanded 
for small 5, and F n was fitted to the resulting Gaussian repulsion between x and y. We would expect 
that the detailed knowledge of the shape and spacing of several individual eigenvalues combined with 
independent measurements at different topologies could provide an alternative way to determine both 
low-energy constants. 



5 Conclusions 

We have derived new analytical expressions for the probability distribution of the fcth lowest Dirac 
operator eigenvalue in QCD with three (or more) colours with spontaneous breaking of chiral sym- 
metry. Extending earlier results we have shown how these distributions become modified when the 
Dirac operator couples to imaginary chemical potential, for arbitrary topology and including the par- 
tially quenched case. Because the deformation parameter of the spectrum is jl = [iF^yfV (where [i is 
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the externally supplied isospin chemical potential), one can use this to extract the magnitude of F n 
through lattice gauge theory simulations. This been already successfully done for the first eigenvalue 
at v = 0, after including first order finite- volume corrections via F we tf. Our results thus offer further 
independent checks for higher topology or higher eigenvalues, given they remain in the e-regime. 

There is a clear sensitivity to F n in these individual eigenvalue distributions, but it requires quark 
masses to be quite small. Because present-day lattice gauge configurations may tend to be available for 
relatively large masses, far from the scale of eigenvalues considered here, the eigenvalue distributions 
will look close to quenched and there is then little variation in the distributions as \x is introduced. 
For this reason it would be extremely helpful if one could derive analogous analytical expressions for 
mixed or conditional individual eigenvalue distributions, say to find simultaneously the first eigenvalue 
of T>\ at x and the first eigenvalue of T>2 at U- In the same way as the mixed spectral two-point density 
correlation function p(x, y) considered previously in the literature, that develops a delta-function in 
the fi — > limit, this quantity should lead to a more dramatic numerical signal as chemical potential 
is turned on and thus to an easier way of measuring F n by means of such a technique. It is an open 
challenge to find new mathematical tricks that will be needed to perform this extension of our results. 
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A Identity for Laguerre Polynomials 

In this Appendix we derive the form of the new polynomials given in eq. (|3.17j) . They first appear in 
the manipulations leading to eq. f)3. 15|) in the following form 

„ (M 2, _ A (-r%!(n+l)! f r l f „(Mj\ .^/A 

= H" n! it, TT~y L ^ (~r) Ll nM~ x k) 

1=0 v ' v y 

j=0 y ' t=0 \J ' /■ ■ 

The first line originates from generating polynomials L\ in the lower rows in eq. f)3. 15|) . Next we have 
swapped the sums, Y^j=o Xw=o ~~ ^ SILoSjMj so that the sum in —X% gives L^(— Then we 
use the following identity: 

L n(^) = E {n -mJ+ j)\ w3{l ~ w)n ~ 3 ^ {z) ■ (A - 2) 

to eliminate the argument 1/t from the first Laguerre polynomial, choosing z = M| and w = 1/t. 
The newly appearing sum is swapped again to ^ "_ Y^i=j ■ The remaining sum can be simplified with 
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the help of the identity 

(-Y(n + u-m + i)\ rn _ m+l+1 

=o 



T-V( x \ - V VI I " ~ I ^ rn-m+i+1/ n /* ^ 



with m = n — j and x = — X^. 

The polynomials q^(x) and their derivation differ slightly from the polynomials Q n (x) in [13j. 
First, we allow for non-zero topology here, in contrast to there. Second, in the related formula to 
([3.15P the Laguerre polynomials l? n have to be generated there, instead of L\ here. 

For completeness we also define the generalization of the polynomials Q n (x) that are needed to 
complete the computations of |13| for non-zero topology: 

Qn{M2) = ^HiW^M 11 \~) { ~ S) \i 

n 1 

= (-) nnl T,7T^ LV j( M 2) L n-j\-S 2 ) , (A.4) 

j=0 { > 

with S 2 = +Ntc\s 2 . The derivation goes along exactly the same lines as above and we skip the details 
here. The single difference compared to eq. (I3.17j) is the shift in the negative index from — v — > —v— 1 
in the final result. In Appendix [C] below we give the corresponding expressions for the first eigenvalue 
from 1 13 1 for non-zero v. 



B Integration Theorem over Determinants 



In this Appendix we compute averages of the product of two determinants of different size that both 
contain different variables and different polynomials. The Theorem we give is a slight generalisation 
of Theorem 1 in [26j . It was already used in [13| but since the form we need here is more general we 
provide a short derivation. In particular, it generalises Dyson's Theorem, see Theorem 5.1.4 in [25J, 
valid only for determinants of the same size and of the same polynomials. 

Suppose we have two sets of (bi)orthonormal functions (this includes the case of polynomials 
orthogonal w.r.t. a given weight, since the weight can be absorbed into the functions) (pj=o,i,...( z ) an d 
<Pj=o,i,...( z ) 

dz tpj(z)<pi(z) = 6ji for j, / = 0, 1, . . . . (B.l) 



Furthermore we assume that we have two matrices (independent of z%, . . . , zn 
then compute the following average over N variables Zj with k > 1: 

^4-1,0 • • • Ai^+k-i Bifl ■ ■ • -B^Tv+fc 



Aij and Bij. We can 



X = 



N 



Y[dzj 



A k ,o 
<Po( z i) 



Ak,N+k-l 

yjv+fc-i(^i) 



<A)0*i) 



<Po{ z N) 



Bk',N+k'-l 
<PN+k'-l( z l) 

<PN+k'-li z N) 



(B.2) 



■ fN+k-lK^N) 

where without loss of generality we assume k > k'. We need to define 4 different objects that we refer 
to as kernels 

N+k'-l N+k'-i 
K N + k'ihj)= Yl A iJ B 3J ' K N+k'(^ Z ) = Yl A i,l<Pl( z ) ' 



1=0 



1=0 
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iV+Jfc'-l AT+jfc'-l 
K I N I l k ,(zJ)= J2 WW%.C(^)= E wWl^)- (B.3) 

Only the last kernel has the self-reproducing properties required for Dyson's Theorem to apply: 

dz'K^+y (z, z')K$^ k , (V, z") = Kj^ +k , (z, z") , 

dz'w{z')K^ +kl (z', z') = N + k' . (B.4) 

Obviously the contraction of mixed kernels leads to mixed kernels, 

dJKtf+xfa^Ktflxtf J) = K I N+k *(i,j) , 
dz'K^ii, z')K^{z', z) = K$ +W {i, z) , 

dz'K l ^ k ,{z,z')Km kl {z\i) = . (B.5) 

The result is now given in terms of the matrix 



By 



-^AT+fc' (hj) i — 1) — j k, J ' — 1, — , y 
Ai,j+N-i i = l,...,k,j = k' + l,k' + 2,...,k 



The claim is that the integral ()B.2[) is given by 

X = (-) N ^- k ^N\det[B] . (B.6) 

By choosing Ay and Bij to be appropriate polynomials of the unintegrated variables, letting (f n (z) 
and <p n {z) be Laguerre polynomials times the weight factor, and choosing N — k variables z\, . . . ,zn — > 
Zk+i, ■ ■ ■ , Z]\f we arrive at the integral of ()3.19p . Taking the integration norm into account, we see 
that (l3~23j) . (|3~24jl and ([3J25D follow from (jR6j) . 

The proof of the formula is by induction on k. We first consider the base case k = k', where 
(IB.6P is just A^! detlK 1 (i, j)]. To derive this, we first transpose one of the matrices in (|B.2p . and then 
multiply them inside a common determinant. The resulting matrix will then consist of blocks with 
each of the four types of kernels. Using the self-reproducing property of K VI and Dyson's Theorem 
|25j . we can carry out the integration, yielding the stated result. 

For the step k — > k + 1 , we expand the first determinant of (|B.2p in the last column: 

k+N+l oo / N \ 

E (-) k+N+1+l / [J d Zj (M A \ k+N+1 det[m A (i)] det[M B ], 
i=i Jo \j=i J 

with Ma {Mb) the first (second) matrix of ()B.2p . and m A («) being Ma with the row i and column 
k + N + 1 removed. The integrands with i > k + l are proportional to (fk+N( z i~k~i), but since there is 
no corresponding (pk+N{ z i-k-i) factor in the expansion of the Mb determinant, these terms are killed 
by orthogonality. Using the induction hypothesis on the remaining terms (proportional to Aj fe+jy) 
and comparing to the expansion of deti3 in the last column, we verify (|B.6p . 
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C Non-Zero Topology for (x): an Alternative Formulation 

In this appendix we derive an alternative formulation of the distribution of the first eigenvalue for 
v > 0. It is based on an earlier paper |13| to where we refer for more details, explicit results were 
given there only for v = 0. Because the calculations in the main part of this paper are quite involved 
it is useful to have an independent, equivalent result as a cross check. 

In [13] one computes first the gap probabilitj@, that the interval [0, s] is empty of all eigenvalues 

Xi . 



i 



7 ' 



dx 1 ...dx N dyi...dy N vl Nf \{x},{y};{m 1 },{m 2 }) . (C.l) 



From this the distribution of the first eigenvalue follows by differentiation, Pi f> (s) 
In steps very similar to the main body of this paper we obtain 

1 



-d.d N '\s). 



E, 



s) 



-N 2 T Cl S 2 



7 K 



(iVi+ATa) 



A; 



N 1 +u(m'/ i )A N2 (m 2 h ) 



dzi 



,dz N e~ N ^^ 



(C.2) 



I2> 



L {M[ 
Lo(M' N \ +u ) 



Ln+N!+u-i(Zi) 



N 2 ' 



L (Z_ 



N) 



Qn+n 2 -i 
Ln+n 2 -i(Zi) 



Ln+n 2 -i(Zn) 



Lq(Zn) ••• Ln+N!+u-i{Zn) 

where we have suppressed all mass independent normalisation factors. The Laguerre polynomials 
Lk(x) = (—) k k\Lk(x) in monic normalisation contain the following masses: 



-NTC 2 m j , j 



1, 



No 



for the flavours of type N 2 , and for flavour Ni we have v additional masses 



M 
M 



12 



12 



-NTdm'j 2 



12 



J 



J 



= l,...,Ni, 
N 1 + l,...,N 1 + v 



(C.3) 



(C.4) 



The degeneracy of the latter can be reinstated by setting ej = at the end of the calculation. 
Alternatively we could write out the resulting derivatives acting on the first determinant in eq. ()C2|) 
for v > 1. Last but not least we have introduced a set of new polynomials containing now generalised 
Laguerre polynomials, as defined in (|A.4p . This completes in principle the computations of |13| for 
non-zero topology. The derivation from the first to the last line in eq. (|A.4p goes along the same lines 
as the previous appendix [A] for the polynomials ^(x), apart from the slightly different identity to be 
used: 



L 



-u-l 
n-j 



(x) 



i=0 



n-j-i 



(X) 



(C.5) 



(j + u)\i\ 

Applying the Theorem derived in the next appendix |D] the integrals over variables Zk can now 
be easily performed. Instead of giving the most general result, distinguishing between different cases 
depending on the numbers N\, N 2 and u, we give a few simple examples for illustration. These can 
be stated in terms of the above polynomials as well as the following kernel: 



/C^(M( 2 ,M 



N-l 

E 

3=0 



L G AM> 2 )Q»{Ml 



(C.6) 



It is denoted there by Eo t o(s,t = 0). 
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Note that only the second polynomial Q V j gets modified when v 7^ 0, compared to [13]. The kernel is 
obviously not symmetric in its arguments. 

In the first example we choose N± = N2 = 1 to get 



4L + o 1} ( s ) 



-N 2 TC V S 2 



12 



-N 2 T Cl S 2 



z[ 1+1 ^ (mi; 777.2) ( m i 



/2 -s 2 ) 



/C^=\(M{ 2 ,M|) /C^=\(5 2 ,M|) 
(M( 2 ) ^ +1 (5 2 ) 



■^jv+i 



(C.7) 



etc., with more rows with Laguerre polynomials for higher v. 

The second example which is probably most relevant for applications is partially quenched, with 
Ni = and N2 = 2. Here the size of the determinant does not grow immediately as the masses 777^2 
of flavour N2 get paired with those generated by topology, v = 1, 2 in our examples: 



E. 



-N 2 TC 1S 2 



z, 



f +2) ( m i,2)( m i - m D 

e -N 2 rc^s 2 



Z 



' +2 '(™i,2)( m i — m|) 



q» n =°(m 2 ) Qr°M) 

)C» N = + \(S 2 ,M 2 ) K" N %(S 2 ,Mi) 



,(0+2) 



(mi^Xm 2 — m 2 )2s 



N+l 



/C^= 2 i(5 2 ,M|) 



(Ci 



In all cases one further differentiation with respect to s yields the distribution of the first eigenvalue. 
C.l The large-TV limit 

We will sketch here first how to take the microscopic large-TV limit of the main building blocks: the 
generalised Laguerre polynomials a , the new polynomials Qj and the kernel /C^. Then we will give 
two explicit examples for the distribution of the first eigenvalue at v > 0. 

The scaling limit is exactly as described in the main body of this paper (and in [13]). so we can 
be brief. We begin with the new polynomials Qj. The sum can be replaced by an integral except for 
a few terms, and we need the following ingredients in the limit N — > 00 together with \i\ 2 — > 



lim (1 — r) 

N,j,k—too 



exp 



where t = j/N, r = k/j , 



lim rci 2 

N-too 



1 . 



(C.9) 



For the Laguerre polynomials the following scaling holds: 



lim L V AM 2 ) 

N,j -too J 



N u (rt) u / 2 (2/m) u I„(mVrt) , 



lim L 

N,j,k— too 



-v-l 
j—k>u 



(S 2 ) 



N~ 



-u-\ 



_ r))-(- +i y 2 (s/2Y +1 I u+1 (sVT^) 



(CIO) 
(C.ll) 



In comparison to v = 0, where we had to split off the s-independent part Lq 1 = 1 from the sum in eq. 
(|A.4p . we now have to separate a total number of v + 1 terms Lq v from the sum. In the large- iV 
limit these terms will seem to be of higher order, but after some cancelations taking place they will 
give a contribution of the same order as the sum itself. 
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To illustrate this we give the first two examples (including the known v = case) for the limiting 



new polynomials lim^oo jT-Q'j (M 2 ) = Q v s (rh; t): 



dr 



o 



-I (mVrt)I 1 (s y /{l-r)t)+ e^ 2 I (Vi 



m) 



(C.12) 



2 /■! 



S 

2m 



1*2 



dr 



rti2 

rt e 2 
1 — r 



Ix{my/H)I 2 {sy/t{T^r)) 



I (mVi) + 



rn 



?2 S 



"I < - 



2m 



1*2 

+e2 d 



dr —I 2 (rhVri)I 3 (sy / t(l - r)) 



,1 — r)2 



Jo(rWi) + ^Ii (Wt) ( ^ + 2<5 2 ) + -^himVt) ( ^ + ^- + <5 4 



S 4 ^2 



The fact that we get the following relation lim Ml 2 ^o( — ) n Qn(^l)/ n ' = -^n( — ^( m 2 + s2 )) fr° m e Q- 
(|A.4j) in the limit of vanishing chemical potentials, implies a series of generalised Sonine identities, 
e.g. for v = 1: 



2^' ' * 



dr-^-/i(mVr)J 2 (sV^ - r) + ii(m)^ + J (m) = I (Vm 2 + s 2 ) 
1 — r J 



(C.13) 



etc. The general identity obtained in this way is stated in the following appendix and proven by 
induction. 

The asymptotic limit of Q u s then easily translates into the limiting form of the microscopic kernel: 

lim i/C^(M?,Mf) = /C£(mi,m 2 ) = \ f dth{m x Vt)Q v s {m 2 -,t) . (C.14) 
iV->oo JS 2 Jq 

We are now ready to give some explicit examples for gap probabilities at non-zero topology. The 
first example is for one flavour, the simplest partially quenched case, with N% = and N 2 = 1 with 
mass rh. While for u = we had 1131 



our new results for topology v = 1, 2 are 

exp 
exp 



1 ~2 1 ?2 

-s z — <r 
4 2 



I (m) 



<ii(*) 



1 ~2 1 ?2 
TS ^ 

4 2 

1 ~2 1 $2 



n? 



£s 1 (s,m) , 



Ii(m) 



shim) 



io(S) 



(C.15) 

(C.16) 
(C.17) 



Note that the normalising partition function is chosen such that it does not vanish at zero mass, 
e.g. zl° +1 \m) = I u (m)/m v '. It also gives rise to an extra factor e^ 2 per unpaired flavour. The 
corresponding first eigenvalues are shown in figure [H together with the corresponding densities from 

D21= 



p(° +1 )(x) = P r n {x) - exp 



Ju(x) f 1 



I v {fh\) 



dTTeh T2p I u (T mi )J u (Tx) , 



(C.18) 
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where the quenched one-matrix model spectral density reads j2j 

Pl uen (x) = % [Ux? ~ J v -i{x)J v+l {xj) 



(C.19) 

It can be seen that our new results for individual eigenvalue distributions nicely follow the corre- 
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Figure 4: The partially quenched case with Aq = and N 2 = 1, at topology v = 1 (left plot) and 
v = 2 (right plot). The first eigenvalue is compared with the density eq. (|C.18|) . for rh = 1.5 at 
5 = 2.5 (middle blue curves) and for m = 1.5 at £ = 0.5 (right red curves). The latter is almost 
indistinguishable from the Nf = 1 flavour one- matrix model result [14|, [TB] (not shown) , to which the 
two-matrix model reduces at 5 — )■ 0. We also display the quenched one-matrix model density from eq. 
()C.19p at the corresponding value of v = 1 and 2 for comparison (left black curves). This curve is 
approached both when m 3> 1 or 6 S> 1. 

sponding known spectral densities almost up to the first local maximum, for the two different values 
of topology and various values of rh and 5 chosen. 



The second set of examples is the gap for two flavours. First we give the case Ni 
topology v = 1 (and 0): 



1 = N 2 at 



E 



(1+1) (§S 



^v=0 l m l ) m 2 



-JCrVi.Tha) 



2/4 



?-(l+l)^ V™'2 

^=1 (mi,m 2 )(m 1 



/C^ =1 (mi,m 2 

/oK) 



/C£ =1 (s,m 2 ) 



with 



mi;m 2 



(it e iHl2 I v (m^t)I u {m 2 Vi) 



(C.20) 



(C.21) 



2{rhim 2 y j 

The corresponding first eigenvalue distributions obtained by differentiating these quantities with re- 
spect to s are plotted in figure [5] against the corresponding spectral density [12J, 

Jq dttJ v (tx)I v (trhi) Jq dtt e^ t2&2 J u (tx)I u (trh 2 ) 



/o dtte^ P I v {trhi)I v (tm 2 ) 



(C.22) 



Here the difference between 1 + 1 flavours (left) and partial quenching (right) becomes visible. Once 
more the individual eigenvalue distributions agree nicely with the corresponding spectral densities. 

The partially quenched two-flavour case, with N\ = and N 2 = 2 reads as follows in the large- ./V 
limit: 



m\m 2 



rhil 2 (mi)li(rh 2 ) - m 2 I 2 (m 2 )Ii(rhi) 



/Cg =1 (s,mi) 

Qv=i(rh l -t = l) Q»= l (m 2 -t = l) 



/Cg =1 (s,m 2 ) 
w=li 



(C.23) 
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Figure 5: The example of two flavours: N\ = N 2 = 1 (left plot), and the partially quenched case with 
N± = and N 2 = 2 (right plot), both at topology v = 1. The first eigenvalue is compared with the 
corresponding two-matrix model densities, for the same values of parameters: rh\ = 1.5, rh 2 = 3.5 at 
5 = 0.5 (right red curves), the same masses at 5 = 2.5 (middle blue curves), and rh\ = 5.5, 777,2 = 3.5 
at 5 = 2.5 (left green curves). The quenched one-matrix model density from eq. (|C19p at v = 1 (top 
black curves) is shown to guide the eye. The curves at smallest 5 = 0.5 (right red curves) almost agree 
for both flavour settings, being very close to the the Nf = 2 flavour one-matrix model result [14J (not 
shown). However, the other two curves differentiate between Nf = 1 + 1 and Nf = + 2 flavours for 
the same parameter values. 



Note again the proper normalisation of the partition function in masses and 5. The corresponding 
density to compare with was derived in [12] where it is given as a determinant (including the corrections 
in the Erratum) 



pl Nf -° +2 \x) = pl uen (x)- exp 



dtte* 52 * 2 J v [xt)I v (rhit)[ - I u (m 2 )(xJ u+ i(x) + 5 z J u (x)) - m 2 4 + i(m 2 )J v (x) 



x (7711/^+1(7711)1^(7712) -7712/1/(7711)^+1(7712' 



(C.24) 







+ J dtte'2 J u {xt)I u {m 2 t)yI u (rhi){xJ u+ i{x) + 5 J v {x)) + m\I v+ i{rhx)J v {x) 

In all cases where we could compare to the approach presented in the main body of this paper the 
curves obtained from the two equivalent approaches agree. 



D Sonine Identities 



In this appendix we prove some integral identities relevant for the appendix O that have to be satisfied 
and that can be reduced to the so-called Sonine integral eq. (|D.1|) below. 

We mentioned already that as a check in the limit of vanishing chemical potential the new polyno- 
mials Q v n have to reduce to the Laguerre polynomials at shifted mass, in order to reproduce the known 
one-matrix models results. While for finite- iV this can be done using the first line of the definition e< 
(IA.4jl . this is not so easy after taking the large- N limit. In fact from eq. (|C.12j) at 5 = and t 



eq. 
10 



The case t 7^ 1 is easily reestablished when rescaling rh,s—¥ yim, y/is. 
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we get the following identities for v = 0, 1, 2: 

To = s f dx X I {mx)I 1 {s^/l-x 2 )+I {m)- I (y / m 2 + s 2 ) = , (D.l) 
Jo vl — x 2 

Xi = — I dx—^h(mx)I 2 (sVl ~ x 2 ) + ^—h{m) + I (m) - I (Vm 2 + s 2 ) = , (D.2) 
m Jo 1 — x z 2m 

f X^ I 5^ S 2 I 

Z 2 = — n / ete 3-I 2 (mx)/ 3 (sVl - x 2 ) + —zh{m) + — h{m) + I (m) - 7 (v m 2 + s 2 ) 

m 2 J (i _ x 2)f 8m 2 2m 

= . (D.3) 

The first equation is the known Sonine identity, see e.g. [27], [28] as well as [13] for an independent 
derivation. It is easier to state and prove the difference between two consecutive integral identities of 
this kind: 

= Zv-Zv-i = — - / dx- — TT I u (mx)I u+ i(sy / l-x 2 )+ - I v (m) - 

m v Jo (\ — x 2 )^~ 2 v v\m v 



s u ,1 



dx- — ir I v -. 1 (mx)I u (sy/l-x 2 ) . (D.4) 

m u 1 Jo (1 - x 2 )2 

Because the way to prove the induction start I\ — Iq = and the induction step is the same, using 
integration by parts, we will be brief. If we use 

*(^) = ^. (■«) 

upon choosing y = s^Jl — x 2 we can write 

xl v+ i{s\/l - x 2 ) ( I u {s\Jl - x 2 ) \ 

(l_ x 2)(,+i)/2 - °* y s(1 _ x 2y/2 J ■ W 

We can therefore rewrite the first integral in eq. (|D.4p as follows: 

1 x u+1 

dx ^-I u (mx)I v+ x(s\J I - x 2 ) = 

o (1 — x 2 ) 2 



x u I v {mx)I v {s^/ 1 — x 2 ) 



s {i-x 2 y/ 2 

I v {m) ( s\ v 1 m f 1 



+ / dx ~^ ^7FT X mlv~i{mx) 

Jo s(l-x l ) v i l 

s\ v 1 m f 1 , , . / r s 

2/ z/! s J (l - x 2)2 



where we have used a Bessel identity pl u (p)' + vl u {p) = pl u -i(p) as well as the series representation 
of the Bessel function to determine the limit at the upper bound x = 1. Inserting this for u = 1 and 
then general ^ yields the induction start and induction step. 
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